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Probing the electron-phonon coupling in ozone-doped graphene by Raman 

spectroscopy 
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We have investigated the effects of ozone treatment on graphene by Raman scattering. Sequential 
ozone short-exposure cycles resulted in increasing the p doping levels as inferred from the blue shift 
of the 2D and G peak frequencies, without introducing significant disorder. The two-phonon 2D 
and 2D' Raman peak intensities show a significant decrease, while, on the contrary, the one-phonon 
G Raman peak intensity remains constant for the whole exposure process. The former reflects the 
dynamics of the photoexcited electrons (holes) and, specifically, the increase of the electron-electron 
scattering rate with doping. From the ratio of 2D to 2D intensities, which remains constant with 
doping, we could extract the ratio of electron-phonon coupling parameters. This ratio is found 
independent on the number of layers up to ten layers. Moreover, the rate of decrease of 2D and 2D' 
intensities with doping was found to slowdown inversely proportional to the number of graphene 
layers, revealing the increase of the electron-electron collision probability. 

PACS numbers: 63.20.kd, 63.22. Re, 78.30.-j, 78.67.Wj 



I. INTRODUCTION 

Graphene linear carrier dispersion in the vicinity of two 
inequivalent points (K, K') of the Brillouin zone creates 
the conditions for the occurrence of unusual effects on the 
dynamics of both electrons (holes) and phonons, which 
are related to the electron-phonon (e-ph) interactionji"— 

In graphene, doping can be tuned by means of the 
field effect, i.e., electric charge induced by capacitive cou- 
pling.— Moreover, being a system entirely exposed to its 
environment, modification of the carrier concentration in 
graphene can follow from molecules adsorbed on the sur- 
face by charge transfer i^ The control of the carrier con- 
centration allows the study of electron-phonon coupling 
(EPC) effects since the e-ph interaction is directly mod- 
ified by changing the Fermi energy level. Raman scat- 
tering measurements in field-effect devices showed the 
dependence of the G peak position and linewidth with 
doping, where G is the one-phonon mode at the T point, 
unveiling tunable optical phonon anomalies i^^— The pos- 
sibility to monitor doping-oriented studies in graphene by 
Raman spectroscopy has provided the basis for a large 
range of applicationiiir— 

Besides the effects of doping on the frequency and 
linewidth of the Raman G peak, the intensity of the 
2D two-phonon signature, or the ratio of peak intensi- 
ties I{2D)/I{G), was found to decrease with increasing 
doping- and used as a tool to qualitatively establish the 
presence of charged impurities J^ii^ More recently, an un- 
derstanding of how the two-phonon Raman peaks inten- 
sity depends on doping has been provided based on fully 



resonant processes 



.17.18 



and their dependence on the 



ph and electron-electron (e-e) collision rates established. 
In Rcfs. 19| and 120. the e-ph scattering rate was not en- 
tirely obtained from experiments as the analysis of the 
2D peak intensity requires the calculation of the e-e scat- 



tering rate. 

Ozone treatment is considered as a promising route 
to enhance the otherwise weak chemical reactivity of 
graphitic structures i2i"— The conductivity of carbon 
nanotubes (CNTs) was found to increase at low ozone 
dose or exposure timei ^^'^^ The proposed mechanism for 
this increase was the ozone adsorption on the CNT sur- 
face, which induces charge transfer effects. At high ozone 
dose (or exposure time), the generation of structural 
modifications and defects seems to be the cause of a re- 
duced conductivity.^^ In this paper we investigate the 
effects of sequential ozone treatment cycles on graphene 
flakes by analyzing the Raman spectrum. Our studies 
show that graphene displays similar changes with ozone 
as those reported in CNTs, i.e., p-type doping without 
introducing significant disorder at low exposure, and dif- 
ferent degrees of bond disruption and surface etching at 
high exposure. Here, we restrict our study to low ex- 
posure conditions and a full account on graphene oxi- 
dation by ozone will be published elsewhere i^ Raman 
spectroscopy tracks the process of gradual p doping of 
the samples, as concluded from both the position and 
the intensity of the Raman peaks. From the latter we 
could determine the EPC for the phonon modes near the 
K point, and to monitor the e-e scattering contribution 
with increasing charge concentration as well as with the 
number of graphene layers. We found a good correlation 
between the rate of the intensity decrease upon doping 
and the number of graphene layers. 



II. SAMPLES AND EXPERIMENTAL METHOD 



Graphene sheets, prepared by micromcchanical cleav- 
age of highly-ordered pyrolytic graphite (HOPG)'', were 
deposited on Si wafers with 300 nm thick thermal silicon 




FIG. 1. AFM images obtained from the same graphene sheet 
before ozone treatment (loop 0) and after five exposure cycles 
to ozone (loop 5). 



for their activation. Thus, these peaks are only detected 
when carbon planes present structural imperfections.— 
Notice the undetectable presence of the latter modes in 
the spectrum of the pristine sample [cf. Fig. [21[a)(loop 
0)] indicating high structural order. 

The peaks denoted as 2D and 2D' in and D' peaks, 
respectively, involving two-phonon processes with oppo- 
site wavectors, which do not require the presence of de- 
fects for their activation. The strongest and featureless 
2D peak in monolayer graphene evolves to a structured 
lineshape as the number of layers increases, revealing the 
electronic band structure, which, in turn, depends on the 
number of stacked layers;^ 



oxide. Ozone treatment of the samples was performed 
at room temperature with a Novascan UV-Ozone clean- 
ing system. The morphology of the graphene sheets was 
studied using a Nanoscope AFM in tapping mode. Micro- 
Raman measurements were carried out at room temper- 
ature in backscattering geometry using a T64000 Jobin- 
Yvon spectrometer with a cooled charge-coupled-device 
detector. In micro-Raman measurements the light was 
focused to a spot with diameter of about 1 /xm through 
a lOOx objective. The 514.5 nm emission line of an Ar+ 
laser was used for excitation with a typical power of only 
120 ^W, in order to prevent structural damage of the 
sample surface by the laser irradiation^^ Raman peak 
lineshape and visible light absorption were used to deter- 
mine the number of layers.^* 

Graphene samples were placed on the ozone cleaner for 
cycles of fixed duration. After each exposure, samples 
were analyzed morphologically by AFM and by Raman 
spectroscopy. The exposure time of 2.5 min during each 
cycle ensures the graphene surface quality is preserved. 
In Fig. [T] AFM images taken from the same sample region 
before and after five exposure cycles [Fig. [IJloop 0) and 
(loop 5), respectively] reveal smooth surfaces. Here, loop 
refers to no exposure to ozone while loop n refers to a 
sample exposed to n consecutive loops. This observation 
is supported by Raman scattering measurements show- 
ing that the structural disorder-related D peak is barely 
detected, which indicates a low amount of damage of the 
graphene surface after being exposed to ozone. 



III. RESULTS AND DISCUSSION 

The Raman spectrum of graphitic structures shows 
common features, i.e., the main one-phonon G peak [cf. 
Fig. IH^a)], and defect-assisted one-phonon D and D' 
peaks, the expected frequencies of which are indicated 
by arrows in Fig. ^a.) .— The G and D' peaks corre- 
spond to phonons at and near the Brillouin T point {Eg 
mode), respectively. The D peak comes from phonons 
near the K points (Ai mode). Both D and D' peaks, 
which are evidence of inter- and intra- valley double res- 
onance processes, respectively, require defect scattering 



A. Second-order Raman peaks intensity 

With increasing ozone exposure, the two-phonon 2D 
and 2D' Raman peak intensities show a significant de- 
crease [see Fig. ^h) and inset to Fig. [2l^b)] while, in con- 
trast, the one-phonon G Raman peak intensity remains 
unchanged with ozone doping. 

The activation of 2D and 2D' peaks involves four-step 
processes where all the states are real, and require en- 
ergy and momentum conservation at every elementary 
step, which means that both two-phonon Raman pro- 
cesses are fully resonant. As a consequence, two-phonon 
Raman spectroscopy is sensitive to the dynamics of the 
photoexcited electron- hole (e-h) pairs, i.e., other than 
e-ph inelastic-scattering processes they can undergo, for 
example, e-e collisions.— Assuming that these two pro- 
cesses are the main scattering mechanisms, the integrated 
intensities over the full linewidth of the 2D and 2D' Ra- 
man peaks, which represent the probabilities of the re- 
spective Raman processes, can be expressed as>^ 



A{2D) = CijKhr 
A{2D') = C"(7r/7)' 



(la) 



(lb) 



where C and C" are constants and 27 denotes the inelastic 
scattering rate of the photoexcited e-h pair written as: 



7 = le~ph +7e-e, 



(2) 



and the phonon emission rate je-ph includes phonons 
near T and K, 



7e-p/i = 7r +7K- 



(3) 



Since the e-e scattering is dependent on carrier den- 
sity, both 2D and 2D' intensities are sensitive to doping 
level. ^^ The e-e scattering rate, 27e_e, was found^'' to 
be proportional to the Fermi energy, Ep, and up to first 
order in Ep it is expressed as: 



le-e^\Ep\f, 



(4) 



where the proportionality coefficient / depends on the 
Coulomb coupling constant. The intensity of the G peak 
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FIG. 2. (a) First- and (b) second-order Raman spectra of 
single-layer graphene taken from the pristine sample (loop 0) 
and after each exposure to ozone up to three cycles (loop 1 to 
loop 3). The inset shows the enlarged spectra in the spectral 
range of the 2D' Raman peak. 



should not depend on doping, since the given Raman pro- 
cess is in off-resonance conditionJS Therefore, the Raman 
peak intensity variation shown in Fig. [51 i.e., decrease 
of two-phonon Raman peaks intensity and unchanged G 
peak intensity, could be attributed to a change of the 
carrier concentration due to ozone treatment. 



B. Frequency shift of Raman bands 

The above interpretation is reinforced by monitoring 
the changes in frequency of the Raman spectrum fea- 
ture after each ozone cycle, given the extensively re- 
ported dependence of the Raman modes frequency on 
doping in graphene. Charge-transfer modification of the 
chemical-bond induces variations of bond lengths (stiff- 
ening/softening), reflected in the variation of the phonon 
frequency. It has been shown that in graphene this effect 
alone does not explain the behavior of the Raman peaks 
frequency with doping, and therefore effects arising from 
the suppression of the Kohn anomaly at T and K points 



must be invoked J. 

The G peak frequency increases for both n and p 
doping, due to the nonadiabatic removal of the Kohn 
anomaly from the F point j^ii Anomalous phonon soft- 
ening, which is seen at low temperature but is smeared 
at 300 K, reflects a resonant e-ph coupling effect when 
the e-h energy gap is smaller than the phonon energy i^ 
Moreover, when the e-h energy gap reaches a value higher 
than the phonon energy a sharp linewidth reduction oc- 
curs as the phonon decay process into e-h pairs suddenly 
ceases /£^ 

The 2D peak shows a different dependence on dop- 
ing, which helps discerning between n- and p-type dop- 
ingi^iiS For electron doping, the 2D peak frequency does 
not change much until a high electron concentration is 
reached, showing a softening for further increase. For 
hole doping, the frequency of the 2D peak increases at 
a rate higher than solely expected from variations of lat- 
tice spacing ..^§J^'^ Although phonon modes contributing 
to 2D peak are away from K points, the effects of the 
Kohn anomaly are not negligible due to the strong e- 
ph coupling. As the doping level increases, the Kohn 
anomaly close to K points is smeared out contributing 
to the stiffening of the 2D peak. In contrast, the influ- 
ence of the Kohn anomaly at T becomes weaker as the 
phonon wavevector departs from this point. Therefore, 
the 2D' peak frequency is expected to be almost insensi- 
tive to changes in doping level. 

Finally, charge concentration changes are not the only 
possible source of phonon spectrum variation. Strain ef- 
fects have already been measured in the Raman spectrum 
of graphene. ^'^ The splitting of the G band into two com- 
ponents displays a shift with applied uniaxial strain at a 
rate of 11 and 32 cm"V%- The 2D and 2D' bands do 
not split and they show a shift amounting to 64 and 35 
cm~-^/%, respectively. 

Figure ini displays the Raman G and 2D peak frequen- 
cies showing a blue shift with increasing exposure to 
ozone. The Raman 2D' peak position, on the contrary, 
displays no variation within the experimental resolution 
(see inset to Fig. [3]) . The larger frequency variation of 
the G peak compared to the 2D peak, together with 
a constant 2D' peak frequency, rules out the effect of 
strain. Therefore, we conclude that the graphene sur- 
face increases its p doping level with each ozone expo- 
sure cycle. The present method to change the amount 
of doping does not require additional processing to fab- 
ricate contacts, which can affect the crystal quality and 
the homogeneity of the properties of the graphene flakes. 
Comparing the measured phonons frequency shift in the 
present study with those found in the literature, we ob- 
tain the overall change of the carrier concentration to be 
Aptot « 5 • 10^^ cm~^. The Raman spectrum after four 
ozone exposure cycles (not displayed) shows no further 
changes, neither in the two-phonon peak intensities nor 
in the peak position of the G and 2D features, indicating 
that the charge concentration, i.e., the ozone adsorption, 
reached a constant value. 
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FIG. 3. Frequency shifts of the Raman 2D vs G peak obtained 
after each ozone exposure cycle. The inset shows the evolution 
of the Raman 2D' vs G peak frequency. 



C. Electron-phonon coupling parameters ratio 

Figure |4] (solid circles) shows the integrated intensity 
of the ID peak, A(2_D) as a function of the integrated 
intensity of the 2D' peak, A(2I?'), both normalized to 
the integrated intensity of the G peak, A{G), calculated 
from spectra of Fig. [21 This dependence, besides showing 
the common decrease of the intensity values, as expected 
from Fig. [2l unveils a constant A(2D)/A(2J9') ratio, as 
they closely follow a straight line. The decrease in in- 
tensity is a direct consequence of the increase of 7e-e 
with doping [cf. Eqs. ([Ta)) and pbl) ]. On the contrary, 
the constant A(2Z))/A(2D') ratio indicates a weaker de- 
pendence of 7e-pft, on doping. This point becomes clear 
taking the ratio of Eqs. (|lap and (|lb[) . which is found to 
be proportional to the square of the emission rates ratio 
of Eg and A\ phonons. 



A{1D)IA{2D') = 2{iKhvY 



(5) 



The linear fit to the single-layer data in Fig. 2] gives a 
value of the slope of 26± 3% and, from Eq. (I5|), we then 
obtain the ratio of phonon emission rates, ^k/it = 3.6. 
This ratio is related to the adimensional EPC parameters 
Ar and Aj^ , as defined in Ref. |l8|, according to 



IK tOout,K Ak 

7r 



^out.r Ar 



(6) 



where uJout,r and ujout,K are the frequecies of the emit- 
ted photons in the respective Raman processes involving 
phonons at T and K point. In order to compare with 
available calculated and experimental values, we relate 
the EPCs parameters to the square of the e-ph inter- 
action matrix elements averaged on the Fermi surface, 
{D^)p and {Dj^)p, as follows^S^ 
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FIG. 4. Experimental A(2D)/A(G) vs A(2D')/A(G') Raman 
peaks area ratio from single-layer graphene (solid circles), 
double-layer graphene (solid squares) and ten-layer graphene 
(solid triangles) spectra measured before ozone treatment 
and after each ozone exposure cycle. The inset shows the 
A(2D')/A(2D) ratio as a function of the number of graphene 
layers, A'^. 



(7) 



where wr and iX)K are the frequencies of the phonons at F 
and K point, respectively. While the value of (_D|j-) of 
graphene and related graphitic structures has been con- 
troversial, due to the lack of reliable experimental data 
of the phonon dispersion around the K point, until re- 
cently j^ the value of (ZJp)^ was obtained^ by relating 
it to the measured dispersion of the Eg mode near Fi^ 
Then taking (-Dr)^ = ^^ {eVlXf, from Eqs. © and 
d?]) we obtain (Dj^) ^ ^ 205 {eV/A)^. The latter is close 
to the EPC value found experimentally from the phonon 
dispersion around the K point of graphite)^ and to the 
computed EPC values of graphene and graphite, when 
nonlocal exchange-correlation effects are includedj5* Ac- 
tually, {D']^) p of graphene and graphite are expected 
to differ slightly due to a large screening effect of the 
exchange interaction in the latter. We measured the 
A{2D)/A{2D') ratio in HOPG and found a value about 
20 % lower than in graphene (cf. inset to Fig. |4]). For 
samples with thicknesses of around twenty graphene lay- 
ers the values of the A{2D) / A{2D') ratio are midway 
between graphene and graphite. The value of {D']^) „ is 
then found to be in the range of 205-188 {eV/XY- 

Finally, concerning the doping dependence of EPC pa- 
rameters in graphene, it has been calculated that (Dj^') 
reduces by « 16% for a variation of Ap = 1.9 • 10^'^ cm~^, 
while (^Dy)p remains unaffected. '^^ Taking into account 
that our estimated Aptot is about four times smaller, the 
expected decrease of (ZJ^^) „ stays within the experimen- 
tal error. 



D. Electron-electron coupling and number of 
graphene layers 

In Fig. m we included data of Raman measurements 
from bilayer and ten-layer graphene deposited on the 
same Si02/Si wafer as the monolayer sample, solid 
squares and triangles, respectively. The experimental 
Raman intensities show a good correlation with the 
monolayer data, as they closely follow the linear rela- 
tion with slope value of 26 (thick line in Fig. |4]). A 
significant feature is the decrease of the intensity rate 
change, given by the difference between intensities in suc- 
cessive treatment loops, which is seen as a higher con- 
centration of data points in Fig. 31 when comparing 
single- layer (SL), bilayer (BL) and ten-layer (TL) sam- 
ples. The dot-dashed line arrows in Fig. 2] are plotted 
to illustrate the data contraction with increasing num- 
ber of layers, which is found to be related by the ratio 
AA(2I?)sL « 2/:^A{2D)bl ~ 10AA(2i:))TL, within an 
error of approximatively 10 %. Since, as discussed pre- 
viously, the decrease of A{2D) and A{2D'), that follows 
ozone treatment, is related to the increase of hole con- 
centration, the result of Fig. 0] suggests that the rate of 
decrease of the intensity, l\A{2D) [l\A{2D')], with dop- 
ing becomes smaller by a factor inversely proportional to 
the number of graphene layers. 

In order to understand the relation /S.A{2D) oc iV~^ 
[/S.A{2D') ex N~^] , we first recall that a Bernal stack with 
N layers, N even [N odd], has N/2 [(iV-f l)/2] electron- 
like and N/2 [{N+l)/2] hole-like subbands almost touch- 
ing at the K point, with a band overlap smaller than 
41 meV. Additional N/2 [{N - l)/2] electron-like and 
N/2 [(TV — l)/2] hole-like outer subbands appear with 
decreasing energy separation with layer number, reach- 
ing a maximum value of 0.4 eV in bilayer graphene. 
The outermost subband is found in the range from 0.4 
eV in bilayer to « 0.8 eV for 20 graphene layers4S As 
the number of layers increases, the number of e-ph pro- 
cesses that contributes to the 2D {2D') integrated in- 
tensity also grows and it is determined by the selection 
rules for optical excitations and for electron scattering 
by phonons.*^ Since the integrated intensities contain the 
weighted probability of the different processes involving 
phonons with close wavevector, the A{2D')/A{2D) ratio 
only depends on the EPC parameters which, up to ten 
layers, remain constant as inferred from the lineal rela- 
tion of the data in Fig. |4l The 7e-e scattering rate, on 
the other hand, experiences an increase as the intersub- 
band e-e collisions are allowed, following the appearance 
of more subbands with increasing number of layers. Al- 
though for the Ep range studied in the present work the 
parabolicity of the energy subbands should be taken into 
account, for illustrative purposes the high doping case is 
considered where the energy dispersion of the subbands 
can be taken as linear . ^'^'^^ Then, in the approximation 
of small momentum transfer, we can generalize Eq. ^ 
to N layers as, 7e-e = N \Ep\ f. Considering hnear dis- 
persive subbands, the carrier concentration, p, is given 



hy p = NuEp. For highly doped samples, the contribu- 
tion from e-ph scattering to the total scattering rate can 
be neglected compared to the much bigger e-e scattering 
component. Thus, the integrated intensity for multilayer 
graphene can be written as: 



A{2D) 
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c 



7| 



N'P^' 



(8) 



and its rate of decrease upon the change of the carrier 
concentration, Ap, is 



AA(2D) 
Ap 



~C 



ll 



Npp 



2»,2 ■ 



(9) 



which accounts for the contraction of the data in Fig. |3] 
with increasing number of layers. 

Based on Eq. [5J the approximate dependence of the ex- 
perimental two-phonon Raman intensity with the inverse 
of N proves that the amount of initial, p, and transferred, 
Ap, charge concentration in the different samples, SLG, 
BLG and TLG, ought to be similar. Taking into account 
that all samples were placed on the same wafer and ex- 
posed to the same treatments, it is reasonable to assume 
that the extent of unintentional or background doping 
and adsorbate coverage was comparable among them. 
The dependence of the two-phonon integrated intensity 
on the number of layers is two- fold, as seen from Eq. ([8]) . 
First, the charge concentration is distributed among the 
subbands, the number of which increase with the number 
of layers, leading to a decrease of the probability of e-e 
collision. Second, with increasing number of subbands 
the number of allowed e-e processes also increases and, 
with it, the probability of e-e collision. Nevertheless, the 
overall effect is a decrease of the two-phonon Raman in- 
tensity because there is the contribution in the scattering 
rate from the simultaneously excited electron and hole, 
as the square in Eq. ([Ta| reflects. 



IV. CONCLUSIONS 

In conclusion, we have shown that the carrier concen- 
tration in graphene gradually increases with sequential 
ozone short-exposure cycles while preserving its crys- 
tallinity The blue-shift of the G and 2D peak frequen- 
cies is evidence oip doping of the samples. In contrast to 
the G peak intensity, which is found to remain constant, 
the 2D and 2D' peak intensities decrease with increas- 
ing number of ozone exposure cycles, i.e., with increasing 
doping. This effect reflects the responsiveness of the two- 
phonon Raman intensity to the dynamics of photocxcited 
e-h pairs and, in particular, to the contribution of the e-e 
scattering. We used this dependence to extract the EPC 
of phonons near the K point and found a close agree- 
ment with previous experimental and theoretical values. 
We also demonstrated an inverse dependence of the rate 
of decrease of the intensities upon doping on the graphene 



number of layers, reflecting the increased probability of 
e-e scattering with increasing number of layers. 
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